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I. INTRODUCTION 

Many solids, like alloys or doped semiconductors, form 
crystals consisting of two or more chemical species. In 
contrast to amorphous solids they possess, at least ap- 
proximately, a regular lattice whose sites are randomly 
occupied by the different components. Understanding 
their electronic structure is an important task which we 
will address here through an approach that allows for a 
comprehensive description of such substitutionally disor- 
dered systems. This approach, which we call the local 
distribution (LD) approach, considers the local density 
of states (LDOS) Pi{uj), which is a quantity of primary 
importance in systems with prominent local interactions 
or scattering, e.g. in the formation of local magnetic mo- 
ments. What makes the LD approach 'non-standard' is 
that it directly deals with the distribution of the LDOS 
in the spirit that Anderson introduced in his pioneering 
workji While the LDOS is directly related to the am- 
plitude of the electron's wave-functions on a certain lat- 
tice site, its distribution captures the fluctuations of the 
wave-functions through the system. 

The LD approach has been originally developed for a 
description of Anderson localization^, and furthermore 
been applied to systems with topological disorderii^ In 
this article we like to demonstrate that it can also de- 
scribe systems with a bimodal disorder distribution, as 
the binary alloymodel. It has been noted several times 
(see e.g. Refs. yloQ), that the physics of this model is 
only partially covered in a mean field description as pro- 
vided by the coherent potential approximation (CPA). 
Instead one has to account for the specific configuration 
of atoms in the vicinity of a lattice siteiSiS To give a bet- 
ter description, a variety of (cluster) extensions to CPA 
has been devised, which explicitly treat correlations on 
finite clusters (see e.g. Ref. or Ref. 0) . We will show 
that the LD approach can serve as a conclusive exten- 
sion of CPA — indeed it contains the CPA as a limit, see 
App. — which implicitly contains these correlations. 

Besides giving very precise results for non-interacting 
systems, an appealing feature of the LD approach is its 
possible application to interacting disordered systems. 
Recently some progress in this direction has been made 
using the LD approach in combination with dynamical 
mean field theoryiA^ We could e.g. show that polaron 
formation in an electron-phonon-coupled system is en- 



hanced in the presence of impurities, leading to polaron 
like defect statesi^^ The Mott transition in a binary alloy 
is another example, where the LD approach might help 
to substantiate present results. "'^'^ 

The outline of this article is as follows: We will shortly 
introduce the binary alloy model and the associated dis- 
tributions, and then apply the LD approach. As a lim- 
iting case we will consider the (quantum) percolation 
model, and finally conclude. The two appendices contain 
the derivation of the LD approach and its application to 
the Anderson localization problem. 

II. MODEL AND DISTRIBUTIONS 

The simplest model for an electron moving in a crystal 
with substitutional disorder is given by the Hamiltonian 

^ = E^»4c, -^E^s ' (1) 

where t denotes the tight-binding hopping integral be- 
tween nearest neighbour sites on a given lattice, and the 
Ei's are random on-site potentials. 

The model is further specified through the probability 
distribution of the random variables e^. We only con- 
sider models where the are identically independently 
distributed random variables with a fixed distribution 
p(ei). The Anderson model of localization, the binary 
alloy model and the quantum percolation model are ex- 
amples. 

For given values of the e^, i.e. for a specific disorder 
configuration ('one specimen'), the local density of states 
(LDOS) 

Piiuj) = -Im Gii{uj)/iT (2) 

(expressed through the retarded Green function Gii(w)) 
has a definite value on each lattice site i. For the or- 
dered system, translational symmetry implies that pi [uj) 
has the same value for every i. In a disordered system, 
Pi{uj) varies with i, so it takes on different values with 
a certain probability. The corresponding probability dis- 
tribution p{pi,uj) captures the fluctuation of the LDOS 
through the system. It is reasonable to assume that for 
an infinite system p{pi,Lij) does not depend on the explicit 
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values of the e^, but only on the probability distribution 
p{ei). We thus assume that the distribution p{pi,uj) is, 
in a slight 'abuse of language', self-averaging: It does not 
depend on the specific disorder configuration looked at, 
but takes on a definite value in the thermodynamic limit. 
Moreover p{pi,uj) has a second meaning. If we look at 
a fixed lattice site i but consider all possible values of 
the ei, the LDOS is a random variable in its own right, 
whose probability distribution is equal to p{pi,uj) as de- 
fined above. So p{pi,u;) gives (i) the probability that the 
LDOS has a certain value on some lattice site, when all 
ei are fixed, and (ii) the probability that the LDOS has 
a certain value on a fixed lattice site i, when the vary. 
Note that p{pi,uj) does not depend on i, in contrast to 
Pi{lj). 

One basic physical quantity that is calculated from 
p{pi,U!) is the arithmetically averaged density of states 
(DOS) 

CXD 

p{uj) = J p^p{pi,uj) dpt (3) 



which counts the number of states at energy uj. For 
the ordered system, when pi^uj) does not depend on i, 
p{pi,uj) = 5{pi — p{io)) is a ^-peak at the DOS. In the 
presence of disorder p{pi,uj) attains a certain width and 
provides additional information on the character of the 
electronic states in the system. A narrow distribution 
corresponds to more or less homogeneous states, when 
electron scattering is weak, while a broad distribution 
refiects strong scattering leading to very inhomogencous 
states. The essential information on the character of 
states at energy w is thus provided by the distribution 
p{pi,uj), which gives a more detailed description of the 
disordered system than the DOS p{uj) alone. 

In particular we can decide whether states at u are lo- 
calized if we employ the precise definition ofp{pi, uj). The 
retarded Green function Gii{uj), hence the LDOS and its 
distribution, is usually calculated for a complex energy 
Lu + IT] with a small positive imaginary part 77, followed 
by analytical continuation to the real axis, i.e. 77 — > 0. 
A finite rj gives a Lorentzian broadening with respect to 
UJ, i.e. a finite energy resolution. For 77 — > 0, the resolu- 
tion increases until peaks in the LDOS (corresponding to 
poles of Gii{uj)) and bands (corresponding to branch cuts 
of Gii{uj)) can be separated. The behaviour of the dis- 
tribution in the limit 77 — *■ is thus different depending 
on the spectral properties of the system. For a spec- 
trum consisting only of discrete peaks — this corresponds 
to localized states — the distribution becomes singular for 
77 — > 0. Contrary to this we get a regular distribution if 
the spectrum is continuous as for extended band ('Bloch') 
states. We use this property in the study of Anderson lo- 
calization (see Add. IB)) . and will use it for the fragmented 
spectra of the quantum percolation model. 

In this article we obtain p{pi,Lu) directly through the 
LD approach (for details see App. 0. It is natural to 
compare the results for p{u!) of the LD approach to the 



corresponding CPA results. While the lattice enters the 
CPA calculations only through the lattice DOS for the 
ordered system the LD approach is explicitly constructed 
for a Bethe lattice^- . Owing to the absence of loops the 
Bethe lattice is a kind of mean-field approximation to 
(hyper-) cubic lattices. Its particular geometry does not 
affect the qualitative behaviour away from the localiza- 
tion transition, which indeed is similar to a cubic lattice 
in three dimensions, see e.g. our discussion of the quan- 
tum percolation model (Sec. |VJ or the phase diagram 
obtained for the Anderson model (App. 01. Conversely 
the existence of different phases and transitions in be- 
tween is correctly described by the LD approach for di- 
mensions > 3, similar to the success of dynamical mean 
field theory for interacting systems in high dimensions. 
Of course at the phase transition, when the correlation 
length diver ges, results for Bethe and cubic lattices dif- 
fer (see Ref. ilif and App. 0). Weak localization in two 
dimensions, on the other hand, is related to interference 
on closed loops and thus not seen on a Bethe lattice. 



III. BINARY ALLOY MODEL 

The binary alloy model has a bimodal disorder distri- 
bution 

p{e^) ^ CAS{e^ - Ea) + {I - CA)S{e, - Eb) , (4) 

corresponding to a crystal randomly composed of A- 
atoms (B-atoms) at energy Ea {Eb)- We set Ea — 
-A/2, Eb A/2, i.e. I\ = Eb - Ea- Different as- 
pects of the involved physics of the binary a lloy mode l 
have been discussed previously (see e.g. Refs. I5l6ll7ll8l9|) . 
The local electron motion strongly depends on the con- 
figuration of atoms at (and in the vicinity of) a certain 
lattice site. This 'cluster effect', showing up as peaks in 
the DOS, is not obtained by the CPA. We will show that 
the LD approach can give a (more) complete picture with 
fair ease. Especially the DOS becomes accessible and will 
turn out to be important for a thorough understanding. 



A. Low A-concentration — small separation energy 

We first study the situation of low concentration ca — 
0.1 of A-atoms (which is below the classical percolation 
threshold, cf. Sec. 0. The A-atoms can be considered as 
the minority species in a bulk system of B-atoms (doping 
a semiconductor is an example). 

For small separation energy A = 0.3 the energy levels 
of the minority A-atoms lie inside the B-band. The DOS 
mainly consists of the B-band centered at Eb — +0.15, 
but shows some additional structures at the lower band 
edge, which are absent in the CPA DOS (cf. Fig. P). 
They originate from the strong fiuctuations in the local 
environment of an atom, which can be clearly seen in the 
distribution of the LDOS (cf. Fig. QJ. Note the various 
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FIG. 1: Binary alloy with A = 0.3, ca = 0.1. The upper 
panel shows the DOS, obtained within the LD approach and 
CPA. The inset displays the conditional DOS for the cluster 
configuration A — BB, the dashed curve has been calculated 
using Eq. with broadening i] — 10~^ (see text). The lower 
panel shows the distribution of the LDOS at lu — 0.6 (curve 
on top). The six other curves show the distribution of the 
conditional LDOS for a specified cluster of three sites. 

peaks of this distribution (in contrast to the Anderson 
model, see App. 0, and its large width displaying the 
pronounced fluctuations of the electron's wave- functions. 

To understand the consequences of different atoms sit- 
uated at a lattice site we can look at the 'conditional' 
LDOS pf^^(a'), subject to the constraint that an A- re- 
spectively B-atom is located at site i. More generally we 
can specify a certain configuration of atoms on a cluster 
of sites centered at i. By considering larger and larger 
clusters every peak in the distribution of the conditional 
LDOS (and in the DOS) can be attributed to a specific 
configuration (cf. Fig. QJ. As the simplest example the 
pronounced peak at lu ~ —0.35 results from a single A- 
atom surrounded by B-atoms (cf. inset Fig. P). Its ap- 
proximate position and form follow from the simple for- 
mula 

(5) 

for the local Green function of one A-impurity embed- 
ded in a B-lattice. The corresponding DOS, using a 



Lorentzian broadening rj = 0.02 in the energy argument 
LU + irj of the Green function to mimic the effect of tun- 
nelling between different A-atoms, fits the conditional 
DOS very well (dashed curve in inset of Fig.QJ. The com- 
plementary situation of a B-atom surrounded by A-atoms 
contributes to the 'hump' in the B-band (cf. Fig. 
which arises from B-atoms neighbouring to A-atoms. 

It should be noted that Eq. ||SJ) gives a (5-peak outside 
the B-band, corresponding to the impurity state at the 
A-atom, only above a critical value of A. This is another 
evidence that the Bethe lattice should be understood as 
an approximation to lattices in dimensions d > 3. 

B. Low A-concentration — split band case 




FIG. 2: DOS for the binary alloy with A = 2.0, ca = 0.1. 
To broaden the peaks of the 'A-states' we add an artificial 
imaginary part r; = 10~^ to the energy argument lo of the 
Green function Gii{u) in Eq. Q. Arrows mark the position 
of three peaks each in the A- and B-band. 

Increasing the separation energy to A = 2.0 (still with 
CA = 0.1) leads to the split band case (see Fig.[2Jl. The 
majority B-atoms still form a band whose DOS shows 
some additional spikes. 

The concentration of the minority A-atoms is below 
the classical percolation threshold hence only finite A- 
clusters exist (cf. Sec. 0. Due to the large energy 
separation the 'A-states' on these clusters are strongly 
damped through scattering on B-atoms. Accordingly the 
'A-states' do not form a band of (extended) states but a 
fragmented set of peaks with varying height and width, 
reminiscent of the percolation model (cf. Sec. IVj. Again 
these peaks emerge from A-clusters embedded in the B- 
lattice. The central peak &t lu ~ —1.031 whose posi- 
tion can be again calculated with Eq. ||SJ) corresponds to 
a single A-atom surrounded by B-atoms. Similarly the 
two side peaks result from two adjacent A-atoms, and so 
forth. The weight of these peaks decreases exponentially 
as for a ./V-atomic cluster. The complementary con- 
figurations, reversing the role of A- and B-atoms, yield 
spikes in the B-band. Since the concentration of B-atoms 
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is large these spikes do not form peaks but merge with 
the B-band. 

Owing to the different strength of scattering we ex- 
pect that all 'A-states' are localized, while the 'B-band' 
remains extended. In accordance with localization in the 
Anderson model (App. 0) we can extract the nature of 
states from the behaviour of the distribution p{pi,oj) for 
Lo shifted by a small imaginary part r/. Indeed we find 
that, for r/ — !■ 0, the distribution is singular in the energy 
range of the 'A-states' but regular for the B-band, which 
proves that 'A-states' are localized and the 'B-band' re- 
mains extended. 



C. Phase diagram 

In Fig. 13 we show the phase diagram of the binary 
alloy for low concentration ca — 0.1. For A below a 
critical value Ac « 0.5 the spectrum consists of a single 
band, which splits into two bands above A^. The value 
of Ac is a little bit larger than obtained within CPA 
j-^CPA 0.45), but significantly smaller than A = 1, 
when two bands with bandwidth W — 1 centered around 
±A/2 would no longer overlap. 

The phase diagram shows the rich physics of the binary 
alloy As discussed in the previous two sections: (i) for- 
mation of many gaps, leading to a strongly fragmented 
minority 'band', (ii) the bands between the gaps show 
additional spikes, (iii) for split bands and low concentra- 
tions, minority states are localized impurity states, while 
majority states are extended band states. These three 
important effects go beyond the CPA, which gives only 
the band splitting. The CPA works reasonably well for 
the majority band, where the local environment of a lat- 
tice site is less important, but it cannot describe the fine 
structure of the minority band in even a crude way. These 
findings are in full agreement with Ref . 0, and we can ad- 
ditionally show the exact form of the DOS which is one 
outcome of the LD approach. So far we have dealt with 
a low concentration ca of the A-species below the classi- 
cal percolation threshold pc = 1/K for the Bethe lattice. 
Then, when only finite A-clusters exist, strong signatures 
and fragmentation of the DOS could be observed. In- 
creasing K with CA fixed, these signatures are partially 
weakened as long as ca < Pc (see Fig. 0J. Above the 
percolation threshold (ca > Pc), infinite A-clusters exist, 
leading to a band which is no longer fragmented. Never- 
theless the DOS still has peaks. For large K (equivalent 
to Pc Ca) the DOS approaches a semicircular form, 
anticipating the CPA result in the limit K = oo (cf. the 
discussion in the appendices). 

For equal concentration ca = cb — 0.5 and K = 3, 
both atom species are in the percolating regime, which 
corresponds not to a doped material but to a stoichiomet- 
ric compound. The phase diagram (Fig.O shows that the 
signatures in the DOS become more pronounced with in- 
creasing A, when the binary alloy model approaches the 
percolation model (see Sec. Ivj) . 




FIG. 3: (Colour online) Phase diagram of the binary alloy, for 
low concentration ca = 0.1 (upper panel), and equal concen- 
tration Ca ~ Cb ~ 0.5 (lower panel). The solid filled curves 
show the DOS for various A. The dashed curves show the 
CPA band edges, and the dotted lines mark ui = ±A/2±W^/2. 



IV. COMBINED ANDERSON AND BINARY 
ALLOY MODEL 

The binary alloy model shows a distinct tendency to- 
wards peak and gap formation. This behaviour is a 
generic feature of systems with a bimodal disorder dis- 
tribution. For example consider the binary alloy model 
with additional on-site box disorder 7, i.e. 

Pie,) = 2AQ^^/2^\e,-EA\) + ^—^eij/2-\e,-EB\) . 
7 7 

(6) 

The CPA suggests that the system is the combination of 
two rescaled Anderson models. This is consistent with 
its phase diagram concerning the position of mobility 
edges. The LD approach additionally shows that struc- 
tures similar to the 'pure' binary alloy appear (cf. Fig.[SJ. 
Evidently the system is not adequately described in terms 
of rescaled Anderson models. For the parameters consid- 
ered 7 is on the same energy scale as the hopping matrix 
element t (here t = i.e. 7 « 0.56t). It is thus 

reasonable to assume that in a real alloy, e.g. a doped 
semiconductor, the peaked structure of the binary alloy 
DOS can be found. Furthermore, if we extract local- 
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FIG. 4: DOS of the impurity band for A = 2.0 and ca = 
0.2, with increasing number of nearest neighbours K of the 
Bethe lattice. For K = 2,4, the concentration ca is below 
the classical percolation threshold Pc = 1/K, for K = 6,12 
above. In contrast to Fig. |21 we did not broaden the peaks 
with an additional rj. 



ization properties of the system by the rj — )■ 0-limit (of. 
App.^, we find the A-band to be entirely localized (see 
Fig. O . Nevertheless 7 is much smaller than the critical 
disorder (7c ~ 3 x bandwidth) for Anderson localization 
of the A-subband (in agreement with Ref . EB) . Therefore 
localization of the A-band is increased due to the strong 
scattering occurring for the minority band in the binary 
alloy. This suggests that impurity band states in a doped 
semiconductor will almost always be localized. 




FIG. 5: The upper panel shows the DOS for the binary alloy 
(A = 2.0, Ca ~ 0.1) with additional on-site disorder 7 — 0.1. 
(see Eq. |^ . The bottom row displays the dependence on rj 
of the LDOS distribution. Lower left panel; In the center of 
the A-band at = — 1 the distribution becomes singular for 
77 — > 0, i.e. states are localized. Lower right panel: For the 
B-band at w = -1-1 the distribution is independent of 77, and 
has a peaked structure known from the binary alloy model 
(cf. Fig. Ell. 



V. QUANTUM PERCOLATION MODEL 

In the limit Eb — > 00, now keeping Ea — fixed, 
the binary alloy model reduces to the percolation model 
where A-sites are embedded in an impenetrable host 
mediumS^. Certainly, in a real alloy with A < cx), the 
electron has the chance to tunnel through the host bar- 
rier. Nevertheless, the percolation model shows very gen- 
eral features which have already shown up in the split 
band case of the binary alloy model (sec above). 

Let us first study the DOS. If the concentration ca 
of A-atoms is below the classical percolation threshold 
Pc only finite clusters exist. The corresponding spec- 
trum is a pure point spectrum which densely fills the 
energy interval [—W/2, W/2] of the tight-binding band 
('Dirac comb'). An TV-site cluster with occurrence proba- 
bility c^(l — Cyi)^+^ contributes peaks at energies spread 
over the full possible range. Consequently the weight of 
a peak varies in contrast to the Anderson model non- 
monotonically with uj. Different statistical properties of 
the finite clusters can be straightforwardly calculated on 
the Bethe lattice. For instance the weight 

00 

= ^C^-l(l-C^)^+lC^ (7) 

Af=l 

^ fl, CA<0.b 

\il-CA)Vcl, >0.5 ^ ' 

of all finite clusters^^ directly follows by help of the gener- 
ating function for the Catalan numbers Cn, which give 
the number of binary trees with N sites. For concen- 
trations above Pc (= 0.5), when wnn < 1, an infinite 
percolating cluster exists. This cluster can support ex- 
tended states which contribute to an continuous spec- 
trum. Since the 'effective' dimension of the percolating 
cluster is smaller than K (especially close to Pc), the 
bandwidth of the resulting band is smaller than W. 

The DOS shows strong signatures (peaks and dips), 
which have already occurred in the binary alloy model 
(see top picture in Fig. EJ. These signatures arise from 
both isolated finite clusters and finite clusters attached to 
the backbone of the infinite percolating cluster. With in- 
creasing concentration their weight reduces, and the sig- 
natures are washed out. At a; = a pronounced (5-peak 
surrounded by a dip in the DOS exists. With increas- 
ing concentration the peak reduces in weight and the dip 
narrows, eventually both merge to a spike in the band. 

The origin of this central peak and the dip can be un- 
derstood on the same level of reasoning as for Eq. 10 • 
If a single atom is attached to the backbone of the per- 
colating cluster the Green function is modified by the 
additional hopping to this atom. The DOS then shows a 
i5-peak at lu — 0, which is the energy of the state located 
at the atom, and a dip around this peak which arises from 
damping of states on the percolating cluster. The same 
argumentation holds for any finite cluster instead of a 
single atom. Since larger clusters have lower probability 
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quantum mechanical electron scatters off all irregulari- 
ties and is possibly localized. This raises the question of 
a quantum percolation threshold Pq above Pc- 
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FIG. 6: (Colour online) DOS p(<^) for the percolation model, 
for various concentrations ca- Upper panel: Full spectrum 
with Lorentzian broadening — 10~^. Lower left panel: 
Central region around ui = without Lorentzian broaden- 
ing. Lower right panel: Region around the second peak at 
LU — t — l/%/32 « 0.177 without Lorentzian broadening. 



of occurrence the central peak is the most pronounced, 
and exists up to the largest concentrations. 

Nevertheless the percolating cluster is not made up of 
its backbone plus one additional finite cluster, but many 
finite clusters attached to it. The Green function G'^{uj) 
of a (half-) infinite chain with one additional atom at- 
tached to each site obeys the recursion relation 



(9) 



which is solved by G'^{ijj) — G°(tJ — ^^/w), where G^{uj) 
is the Green function of the Bethe lattice. Of course 
G"^(w = 0) = as before, but the diverging real part of 
I/lu produces not a dip but a gap in the DOS around 
uj = 0. The validity of this simple argumentation can 
be tested within the LD approach since its energy res- 
olution is not limited by finite size effects. Without 
Lorentzian broadening, i.e. 77 — > 0, the (5-peaks arising 
from finite cluster states do not contribute to the DOS 
(cf. Sec.^, and only the continuous spectrum from ex- 
tended band states survives. This spectrum indeed shows 
a gap around lu = for sufficiently small concentrations 
(cf. Fig. bottom left). Lowering the concentration, 
gaps open at the energies corresponding to any finite clus- 
ter eigenstate. The first additional gaps open at uj = ±t 
(Fig. El bottom right), corresponding to two site clusters, 
which have the second largest weight among all finite 
clusters. These gaps are filled by peaks from the dense 
spectrum of finite cluster states, thus are absent in the 
DOS for any finite energy resolution (cf. Ref . IT^ . The 
formation of gaps — which is a kind of 'level repulsion' — is 
a significant quantum feature with no counterpart in the 
classical model. 

Besides the DOS the nature of states is important. For 
concentrations ca < Pc, when only finite clusters exist, 
all states are localized and no electron transport is possi- 
ble. For concentrations above Pc a classical electron can 
traverse the system along the percolating cluster. But a 
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FIG. 7: LDOS distribution for the percolation model at lu — 
0.1. Left panel: For ca = 0.7 and two values of regularization 
77 — 10~* and r/ = 10"^". For t] ~* only extended states 
contribute to p{pi,Lu) at finite pi. Right panel: Integrated 
probability p{pi > £,,uj) in dependence on r], for various ca 
and one particular value ^ — 10^'^''. For ca < 0.54 no states 
exist a.t UJ = 0.1, and p{pi,uj) — 0. 

Close to Pc states on the percolating cluster have only 
exponentially small weight, ft is thus difficult to extract 
the information about those states from the background 
of finite cluster contributions. However, we can gain some 
information from the distribution of the LDOS, if ca and 
77 are varied. Again, if we let 77 — > 0, the energy reso- 
lution increases until contributions from extended states 
on the percolating cluster are separated from the discrete 
peaks of localized states (see left panel of Fig. [Tj). For 
fixed energy lo and suitably large concentration ca the 
distribution p{pi,uj) has finite weight at finite pi in the 
limit 77 ^ (see right panel of Fig. [7)), which indicates 
that extended states exist. Their weight can be esti- 
mated through the integrated probability p{pi > £^,uj) — 
p(pi,uj) dpi, with rj < S^. The weight decreases with 
Ca, and abruptly drops to zero at a certain concentra- 
tion (for CA ~ 0.54 a.i ijj — 0.1). Then no states exist 
at (jj, and Pavc(w) — 0. This behaviour has to be at- 
tributed to the fragmentation of the spectrum discussed 
before. At a given energy uj a gap will open for concen- 
trations CA sufficiently close to Pc- Before a gap opens 
at Lo, states are extended. After the gap has opened, no 
states at uj exist. Between the gaps, which open at differ- 
ent concentrations for different energies, extended states 
can survive even for very small ca , although all states are 
extremely damped and practically localized. A definite 
localization transition from extended to localized states 
does however not take place, and the only transition oc- 
curs when the spectrum is fully fragmented, that is at 
the classical transition ca — Pc- We conclude, on the 
basis of these arguments and our numerical results, that, 
for the percolation problem on the Bethe lattice, a quan- 



7 



turn percolation threshold above the classical one does 
not exist. Note that scattering on the percolating cluster 
is of a different type than for the Anderson model. The 
finite clusters attached to the percolating backbone do 
not act as coherent but incoherent scatterers. So states 
on the backbone will not be localized for small ca (i.e. 
strong scattering), or even immediately localized, as in 
one dimension. 



VI. CONCLUSIONS 

In this article we demonstrated how the LD approach 
can be used to study localization and percolative ef- 
fects in alloys. With very moderate computational de- 
mands this scheme suffices to resolve the rich structures 
in the DOS originating from comparably strong disorder 
fluctuations. Even for the extreme limit of the binary 
alloy model, the percolation model, convincing results 
are easily obtained. For instance the question whether 
gaps form in the percolation model could be definitely 
answered, and the possibility of a quantum percolation 
threshold above the classical one could be almost defi- 
nitely ruled out for the Bethe lattice. 

We conclude that the LD approach is a convenient 
framework for investigations of disorder and localization, 
and suggest its application to interacting disordered sys- 
tems. 



APPENDIX A: THE LD APPROACH 

The LD approach has been constructed by Abou- 
Chacra, Anderson, and Thouless'^ on a Bethe lattice. 
There the local Green function G'm(w) = {i\{u! + ir; — 
H)~^\i) can be expressed through Green functions on 
the K neighbouring lattice sites j = 1, . . . ,K, 



Gii{uj) 



K 



(Al) 



where each Gjj{ui) on the r.h.s. of this equation is evalu- 
ated for the lattice with site i removed. Iterating this 
expansion an infinite hierarchy of equations is gener- 
ated ('renormalized perturbation expansion'). -'^^ Instead 
of solving this hierarchy for many particular realizations 
of the Ei and constructing the distribution of Gu (w) af- 
terwards the LD approach manages to solve Eq. (|A1|I 
directly for the distribution. This solution relies on two 
properties arising from the special geometry of the Bethe 
lattice. First, all Green functions in Eq. lAip corre- 
spond to the same geometric situation (one lattice site 
with K neighbouring sites). Hence their distribution is 
identical (although their concrete values differ). More- 
over, removing site i from the lattice, the lattice sites 
j = 1,...^K are unconnected, and the Gjj{Lj) are in- 
dependently distributed. Owing to these two properties 



Eq. IHI) expresses one random variable through K in- 
dependently distributed random variables with the same 
distribution, and can therefore be interpreted as a self- 
consistency equation for the distribution of Gii{uj)mii^ 
Be aware that then the indices «, j do not denote spe- 
cific lattice sites but certain realizations of the random 
variable Gii{uj). 

For the ordered system (e^ — 0) all Green functions in 
Eq. (jAl|) are identical, and the Green function G"(aj) = 
i8/W^){Lo - _ of the Bethe latticeii, with 

bandwidth W = 4t\/K, is obtained. For a disordered 
system the solution of Eq. ljAl|l is obtained through a 
Monte-Carlo procedure (Gibbs sampling). The distribu- 
tion is represented through a sample of typically 10^ up 
to 10'' elements, depending on the respective case stud- 
ied. At each step of iteration a new sample is constructed 
whose elements are calculated through Eq. IjAip with a 
randomly chosen and K elements drawn from the pre- 
vious sample. A hundred up to some thousand iterations 
are necessary to guarantee convergence. The resulting 
computation time on a standard desktop PC ranges from 
few minutes to some hours. 

The LD approach comprises the coherent potential 
approximation (CPA) in the limit of infinite coordina- 
tion number K = oo. Taking into account the scaling 
t oc 1/\/K of the hopping matrix element, the sum over 
J in Eq. (|A1|) can — assuming the central limit theorem 
to be applicable — be replaced by the arithmetic average 
of Gjj, and the CPA is recovered. 

In course of its construction the LD approach works on 
an infinite lattice (no boundaries, no finite size effects, 
no finite energy resolution). In the numerical solution 
the size of the sample determines the resolution which 
the distribution is sampled with. The resolution can be 
easily enhanced by increasing the sample size (see below) . 



APPENDIX B: ANDERSON MODEL 

We describe in this appendix the application of the 
LD approach to the Anderson localization problem. Al- 
though this is not the main objective of this article, it 
will help to underline the generality of the LD approach. 

The Anderson model — which is the prototype model 
showing a localization transition — is given by Eq. ^ for 
a uniform distribution of in the interval [—7/2,7/2], 



I2 - 



(Bl) 



The characteristics of localization show up in p{pi,uj) 
(cf. Fig. (|SJ). Impurity scattering causes the distribution 
to be strongly asymmetric and broad even for extended 
states. As one consequence the DOS p{uj) is on a different 
scale than 'typical' values of the distribution, e.g. the 
most probable value. On approaching the localization 
transition the asymmetry further increases. Much weight 
is transferred to large values of pi, while the 'typical' 
values tends to zero. 



8 



>-2 



1 ' 


o—o LD 






X LD (25M) 






-- TMT 






simpleTMT 


20 




Q. 




X 






\) X 















10= 








\ 




Q. 


\ 






1 , ]'■■■. 


/ 1 , 1 







' mpv 1 


1 

















10* 


10-^ 


10° 








— n=io"° 










11=10° 










n=io"'° 














10 


10 


10 = 



0.2 0.4 0.6 0.8 1 
(0 



FIG. 8: Left panel: Phase diagram of the Anderson modeL 
The soUd (dashed, dotted) curve shows the mobihty edge tra- 
jectory calculated within the LD approach (TMT, simplified 
TMT) by means of the — > limit, using a sample with 5x 10'* 
elements. The crosses indicate points in the (a;,7)-plane cor- 
responding to the mobility edge position for an extremely 
large sample of 2.5 x 10^ elements. Upper right panel: prob- 
ability distribution of LDOS for extended states, at 7 = 1.5, 
io — 0.0. The distribution is - for not too large values of rj 
- independent of rj, but the most probable value ('mpv') is 
some orders of magnitude smaller than the (geometric) DOS 
p{lj) (pgeo('^)). Lower right panel: probability distribution 
p{pi, uj) for localized states, at 7 = 1.5, u) = 0.9, and three 
values of rj. 



A closer look at the distribution reveals that near the 
localization transition it shows a power law behaviour 
over a wide range of pi, with an exponent ~ 1.45, which 
is reasonably close to analytical results obtained from 
field theoretical considerationsii^ 

Passing the localization transition the distribution be- 
comes singular, corresponding to a transition from con- 
tinuous to discrete spectrum. This characteristic change 
reflects itself most clearly in the dependence of the distri- 
bution on the 77-regularization, when Gu^lu) is calculated 
for a complex energy argument Lo+irj. While the distribu- 
tion for extended states is stable when decreasing 77, the 
relevant scale for localized states is entirely set through 
77, and the distribution becomes singular for 77 ^ (see 
Fig. ^ . We can use this different behaviour as a localiza- 
tion criterion, if we perform the 77 — > limit numerically. 
This criterion does not depend on any a priori choice or 
approximation, hence should be considered 'numerically 
exact'. As we have already mentioned (cf. Add. the 
finite size of the Monte-Carlo sample sets the resolution 
for sampling p{pi,Lu). If the sample is too small (and 
thus the resolution too low) an almost singular distri- 
bution will be falsely detected as singular while a larger 
sample correctly gives a regular distribution correspond- 
ing to extended states. Accordingly the mobility edge 
is shifted to larger values of disorder on increasing the 
sample size (Fig.|S}. However, the points on the mobility 
edge trajectory readily stabilize if the sample is chosen 



large enough, and a precise determination of the mobility 
edge is possible. 

The LD approach phase diagram for the Anderson 
model which is then obtained shows the characteris- 
tic features of the localization problem in three dimen- 
sionsiiSiiS These can be most simply understood to arise 
from the interplay of two competing effects: While, for 
small disorder, tunnelling between shallow impurities 
produces extended states outside the tight-binding band, 
strong scattering on deep impurities takes place with in- 
creasing disorder, starting to localize formerly extended 
states. Therefore a reentrant behaviour of the mobility 
edge trajectory and the existence of a critical disorder jc 
(7c w 3.0 ior K = 2 neighbour sites of the Bethe lattice) 
for complete localization of all states is found. 

According to the mean-field type of approximation, the 
critical disorder 7c for the Bethe lattice is larger than for 
a cubic lattice. For K — > 00, 7c grows without bound, 
i.e. localization is absent in ii" = 00 (where the LD ap- 
proach reduces to the CPA) . For K = 1 the Bethe lattice 
is a one-dimensional chain, where all states are known to 
be localized for arbitrary disorder. The sum in Eq. (|A1|) 
then contains only one Gjj (w), which implies that the LD 
sampling scheme cannot converge to a stable distribu- 
tion, since different elements of the sample never become 
related to each other during the sampling. This instabil- 
ity expresses the particular one-dimensional localization 
behaviour. 

For comparison to the LD approach we show in Fig. |S1 
the mobility edge trajectory obtained within a recently 
proposed mean-field like approach to Anderson localiza- 
tion, the so-called typical medium theory (TMT)i^2iSi 
This TMT modifies the CPA by reformulating its self- 
consistency condition in terms of the geometrically aver- 
aged DOS 



Pgoo(t^) = exp / {\npi)p{pi,Lu)dpi 



(B2) 



This average is known to be critical at the localization 
transition, since it puts much weight at low values of pi. 
It nevertheless does not approximate the 'typical' values, 
see Fig. |H1 

If we compare the lines of vanishing Pgco 

{oj) from TMT 

(i.e. the 'TMT mobility edges') with the LD approach 
phase diagram we see the consequences of this modifi- 
cation: (i) the critical disorder predicted is significantly 
smaller, and (ii) the reentrant behaviour of the mobility 
edge is entirely missed. 

Remember that the CPA is obtained in the well defined 
limit of infinite coordination number K — 00, therefore 
is a controlled approximation. In the TMT construction 
some ambiguity enters in the choice of the average used. 
In fact TMT can be further simplified, replacing pgcoi^) 
by Psimp = min{pi(ej = j/2),pi{ei = -7/2)}- This 'av- 
eraging procedure' drastically overestimates the strength 
of impurity scattering, and the 'simplified TMT' is surely 
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far away from any reasonable description of the underly- 
ing physics. However, the phase diagram obtained is very 
similar to the TMT one. Of course the critical disorder 
has to be even smaller than in TMT. 

Apparently TMT captures strong impurity scattering 
which is partially neglected in the CPA. However, Ander- 



son localization is not merely a result of strong scatter- 
ing but quantum interference due to coherent scattering. 
The similarity between the TMT and 'simplified TMT' 
results indicates that TMT might not adequately include 
these interference effects. 
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The (half-infinite) Bethe lattice is a loop-free tree with 
a semicircular DOS p{uj) = {A/-kW'^)^W^ - ico^ii^ If not 
stated otherwise, we use a Bethe lattice with K = 2 nearest 
neighbours to any lattice site. The hopping matrix element 
t is chosen to give a full bandwidth W = 4ty/K = 1, so 
energies will be measured in units of W. 
On the Bethe lattice site and bond percolation are equiv- 
alent. 

Again on a Bethe lattice with coordination number K — 2, 
cf. first footnote. 



